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We have in mind a register of qnbits for an quantum information system, and consider its de¬ 
coherence in an idealized but typical situation. Spontaneous decay and other couplings to the far 
environment considered as the world outside the quantum apparatus will be neglected, while cou¬ 
plings to quantum states within the apparatus, i.e. to a near environment are assumed to dominate. 

Thus the central system couples to the near environment which in turn couples to a far environment. 
Considering that the dynamics in the near environment is not sufficiently well known or controllable, 
we shall use random matrix methods to obtain analytic results. We consider a simplified situation 
where the central system suffers weak dephasing from the near environment, which in turn is cou¬ 
pled randomly to the far environment. We find the anti-intuitive result that increasing the coupling 
between near and far environment actually protects the central qubit. 
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In many quantum optics experiments and quantum in¬ 
formation devices we find the following situation: The 
central system, well protected from simple decoherence 
processes such as spontaneous emission or direct coupling 
to a structureless heat bath, still suffers some decoher¬ 
ence from the coupling to the quantum part of the ap¬ 
paratus. We will call the former the far environment 
and the latter the near environment. In such a tripartite 
system without direct coupling between central system 
and far environment, we find that increasing the cou¬ 
pling of the near to the far environment can protect the 
central system against decoherence. In the setting of 
the Haroche experiment [1, 2] the late M. C. Nemes dis¬ 
cussed this somewhat anti-intuitive fact with one of the 
authors [3] twelve years ago. More recently, additional 
numerical evidence has appeared [4] , some of which were 
master thesis related to the present work [5-7]. Finally, 
it was shown in Ref. [8] that a protected subspace can 
appear in a strong coupling limit. 

We now wish to construct a model that allows some 
analytical treatment and simultaneously has some claim 
to universality. Indeed the intermediate environment in 
a quantum information system typically consists of quan¬ 
tum states that are not used but are unavoidably present. 
While ordered and systematic couplings can be mini¬ 
mized with good technology, it is plausible that the un¬ 
controlled remanent will have a random structure. Con¬ 
sidering its minimum information character [9], we use 
random matrix theory (RMT) of decoherence [10-12]. 
There are several examples, which involve chaotic or ir¬ 
regular dynamics explicitly: e.g. the coupling of two-level 
atoms to quantum systems with classical chaotic ana¬ 
log [13, 14], and the experimental realization of quantum 
chaos in a chain of three level atoms in [15]. The recent 
advances in the design and control of chains of individual 
ions [16], may lead to similar experimental systems. We 
simplify the picture by limiting the coupling between cen¬ 


tral system and near environment to dephasing [17, 18], 
and assume the coupling between near and far environ¬ 
ment to be of a tensor product form; see Eq. (1), below. 
We can then take advantage of analytic expressions that 
exist for dephasing in absence of the far environment [17- 
21] and treat the effect of the far environment within a 
linear response calculation. We thus obtain analytic ex¬ 
pressions for weak couplings between near and far en¬ 
vironment, note that this range of coupling strengths is 
exactly opposite to that treated in Ref. [8] . To study the 
effect of the far environment beyond the linear response 
approximation, we perform numerical simulations, using 
a modified Caldeira-Leggett master equation [22] , whose 
equivalence to RMT models has first been discussed in 
Ref. [23]. 

Model: The full system consists of three parts, the cen¬ 
tral system, the near environment and the far environ¬ 
ment with Hilbert spaces T-Lc, T~Le and 77f, respectively. 
The unitary evolution of the entire system is given by 
the Hamiltonian 

Htot = i7o + ® K ® If + 7 lc 0 K' 0 Vf (1) 

where Hq = /ic 0 le,f + lc 0 Tie 0 If + lc,e 0 Hf- Tracing 
out both environments leads to the reduced dynamics of 
the central system gdt) = tre,f[ptot(i)], with 

gtot{t) = exp{-iHtott/h) Pc 0 6'e,f exp{iHtott/h) , (2) 

where Pc and pe,f represent the initial states of the cen¬ 
tral system (typically assumed to be pure), and the en¬ 
vironment (near and far environment), respectively. The 
couplings are given by the tensor products Uc 0 14 (be¬ 
tween central system and near environment) and 7 
(between near and far environment). The former will be 
chosen as dephasing, such that [hc,Vc\ = 0. Such cou¬ 
plings are frequently used, as they simplify calculations 
and maintain many essential properties. 
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Dynamics: We write the Hamiltonian as Htot = 

= (^6j 14 ) 0lLf H“'7 V^0Vf , (3) 

where the set of states { |j) }j is a common eigenba- 
sis of he and Vc, while Sj and lyj are the corresponding 
eigenvalues. The evolution of the whole system can be 
written as ptot(0 = Y.jkP 3 k{ 0 ) \j){k\ ® where 

Qc = ^jkPjk{0) |j)(^| is the initial state of the central 
system, and 

g0.fc)(^) = gxp ^ pe.f exp (ii?ij t/7i) ■ (4) 

We hnd for the matrix elements of the reduced state of 
the central system: Pjk{t) = pjfc(O) tre,f[p^^’^^(t)]. Since 
trey[pi'^’'’4 (i)] = the diagonal elements are constant in 
time, while the off-diagonal ones (i.e. the coherences) are 
given as expectation values of generalized echo operators 
in the composite environment; see Eq. (6) and Ref. [24]. 

In other words focussing on an individual matrix ele¬ 
ment, Pjkit), we may introduce 

Hx =Ho + XVes = He + y, K , i^o = + ^k K , (5) 

such that A = {vj — Vk) 14. This allows to connect the 
coherences for vanishing coupling to the far environment 
(7 —>■ 0), with fidelity amplitudes [17, 18]. Introducing 
the relative coherences 

f (f\ _ _ Pjk (t) _ 

= tre.f[e-‘^^-‘/'‘pe.fe'^“-‘/'‘] , ( 6 ) 

where Hx,^ = Hx 0 If -I- le ® -f 7 14 ' 0 14 , we find 
that /A,o(t) = f\{t) with 

fxit) = tre[MA(t)trf(pe.f)] , Mxit) = 

( 7 ) 

Hence, fx,o(t) becomes the fidelity amplitude for perturb¬ 
ing the Hamiltonian Hq by A14ff, given the initial state 
trf(pe 7 ) in the near environment. 

Modeling the effect of the far environment, we con¬ 
sider the simplest possible situation, where random ma¬ 
trix and master equation descriptions are equivalent [23] . 
This will allow to use a master equation for numerics 
and the random matrix model for analytics. In Ref. [11], 
it has been shown that without central system, the co¬ 
herences in the near environment decay with the rate 
r = 27r iVe 7 ^/(?i df), which is just A4 times the Fermi- 
golden-rule rate for transitions between individual states. 
Here, 7 ^ replaces the magnitude squared of the coupling 
matrix elements, since in our case 14 is chosen from a ap¬ 
propriately normalized random matrix ensemble, while 
df denotes the average level spacing (or inverse level den¬ 
sity) of idf, and Nf, is the dimension of the near environ¬ 
ment. Choosing the Caldeira-Leggett master equation to 
describe the far environment, one obtains practically the 


same reduced dynamics. The only difference is that T 
then depends on the dissipation constant and the tem¬ 
perature [5, 6 ]. 

Linear response calculation: Applying the linear re¬ 
sponse approximation in the Fermi golden rule regime 
to the coupling between near and far environment (see 
appendix), one arrives at 

fxAt) (1 - rt) fx{t) +T f At fxir) fx[t - r) . ( 8 ) 

do 

The symbol ~ means equal up to O(r^), and from know 
on, we replace 7 by the physically more meaningful de¬ 
coherence rate T. As we will see below, Eq. ( 8 ) equation 
is valid as long as Ft <C 1. It constitutes our main re¬ 
sult. The details of its derivation can be found in the 
appendix. Note that the result is valid for any Hamil¬ 
tonian He, which may result in very different behaviors 
of fx{t)- The only necessary assumption is that I 4 ' is 
sufficiently random. 

Generally, we find that increasing the coupling 
strength to the far environment is indeed slowing down 
the decoherence in the central system. Depending on the 
interaction strength between central system and near en¬ 
vironment, and on the functional form of fxi^), the effect 
can be more or less pronounced. This can be demon¬ 
strated for generic systems, where one often finds that 
fx(t) changes from an exponential decay in the Fermi 
golden rule regime to a Gaussian decay in the pertur¬ 
bative regime [19, 25]. In the former the effect is zero, 
which can also be understood in physical terms. In that 
regime the temporal correlations (I 4 ff(t)Kff(t 4 ) of the 
perturbation in the interaction picture decay very fast - 
on a time scale tcorr ^ tdec(ce), the decoherence time in 
the central system. Therefore, even if the decoherence 
time in the near environment tdec(ef) = F”^ is smaller 
than tdec(ce), as long as tdec(ef) > Aorr, the far envi¬ 
ronment will have no effect on the decoherence in the 
central system. In the perturbative regime by contrast, 
fx(t) = such that fx,r{t) ~ gr/\{Xt) with 

Paix) = {I —ax) e~^ +a y'7r/2 e“^ erf(x/-\/^ . (9) 

Although exact analytical results for fx{t) exist [20, 21], 
for simplicity, we will compare our results to the exponen¬ 
tiated linear response (ELR) expression from Ref. [19]. 

/fL^(t) = exp[-A 2 C'(t)] ( 10 ) 

et' 

C{t) = + TTt — 47r^ / At' / dt" 62 (t") , 

Jo Jo 

where 62 (t) is the two-point form factor [26]. 

Caldeira-Leggett master equation: For full-fledged nu¬ 
merical random matrix calculations, we would need to 
work in the Hilbert space of near and far environment. 
For the far environment, we would need a smaller mean 
level spacing in combination with a larger spectral span, 
as compared to the near environment. Still, in order to 
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FIG. 1. Relative coherence f\.r{t) in the Caldeira-Leggett 
model in the perturbative regime (A = 0.02), for different 
values of the coupling to the heat bath: F/A = 0 (squares), 1.0 
(shaped crosses), 5.0 (inverted triangles) and 10.0 (crosses). 
The time is in units of h/do where 2711%/do is the Heisenberg 
time and do the mean level spacing in the near environment. 
In all cases = 50 and Urun = 1 000 realizations. 


justify the use of RMT, we would need as many levels 
as possible also in the near environment. Such random 
matrix calculations are not viable, due to the dimension 
of the Hamiltonian matrices involved. 

We will therefore use an approach which allows to work 
in the Hilbert space of the near environment alone, tak¬ 
ing the effect of the far environment into account via the 
Caldeira-Leggett master equation [22], where we replace 
the diagonal matrix representation of the harmonic os¬ 
cillator Hamiltonian with a random matrix, defined as 
in Eq. (5). We choose both matrices Hq and Ves from 
the Gaussian orthogonal ensemble (GOE). We scale Hq 
in such a way that the mean level spacing becomes one 
in the center of the spectrum. The matrix elements of 
Veff are chosen to have the variances = I -I- Sij. 

In that way, the strength of the perturbation (implied by 
the dephasing coupling to the central system), measured 
in units of the mean level spacing do, is given by A. In 
the following figures we scale time by the Heisenberg time 
tn = 2t: fi/do- 

Numerical simulations: The use of random matrices 
requires a Monte Carlo average over many realizations. 
As a sufficiently large but still numerically manageable 
dimension of the environment, we choose Ng = 50, and 
perform averages over Urun = 1 000 realisations. The gen¬ 
eral behaviour of the relative coherence f\.r{t) is shown 
in Eig. I. Here, we choose A = 0.02 for the dephasing 
coupling, and different values for the coupling strength 
r between both environments. The figure clearly shows 
that the coherence decays slower by increasing T. 

In the remaining figures, we evaluate the quality of 
our analytical result from Eq. (8). Eor a better quanti¬ 
tative comparison, we subtract the ELR approximation 



FIG. 2. (a) Relative coherence subtracted by the ex¬ 
ponentiated linear response approximation to the fi¬ 

delity amplitude, for A = 0.1 (cross-over regime) and different 
values for the coupling to the heat bath: F/A = 0 (squares), 
0.1 (circles), 0.5 (triangles), and 1.0 (shaped crosses). The 
thick solid lines, show the analytical result according to 
Eq. (8), and the nearest thin dotted lines show three statis¬ 
tically independent ensemble averages for each case, (b) The 
same quantity as in panel (a), but for the theoretical curves 
the values for F are replaced by best fit values Fflt such that: 
Fflt/A = 0.097 (circles), 0.44 (triangles), and 0.77 (shaped 
crosses), obtained for the region 0 < t < 15. The units of 
time are the same as in Fig 1. 

for pure fidelity decay, from both, the numerical 
simulation and the analytical approximation for /A,r(t)- 
Note that for the function f\{t) appearing in the ana¬ 
lytical expression, we use numerical results with much 
improved accuracy. These are obtained from numeri¬ 
cal simulations without far environment and some subse¬ 
quent spline-fitting for facilitating the evaluation of the 
integral in Eq. (8). In that way, the accuracy can be 
greatly improved. The numerical result for f\,o{t) differs 
from the ELR result due to the fact that ff^^{t) is only 
an approximation, but also because of the level density 
over the spectral range of Hq. The remaining difference 
to the exact analytical expression found by Stockmann 
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FIG. 3. Comparison between Caldeira-Leggett simulations 
and the linear response theory as in Fig. 2, but for A = 0.02. 
The units of time are the same as in Fig 1. For the the¬ 
oretical curves, we rescaled Ffit as explained in the text. 
The different cases shown are: (squares) F = 0, circles 
r(rflt) = 0.002 (0.00195), triangles 0.004 (0.0039), shaped 
crosses 0.01 (0.00858), diamonds 0.02 (0.0159), inverted tri¬ 
angles 0.04 (0.0263), crosses 0.1 (0.0457). 


and Schafer [20, 21], is due to the fact that the trace in 
Eq. (7) includes the full spectral range where the level 
density varies according to the semi-circle law. 

In Fig. 2(a) we consider the case A = 0.1, which is in 
the cross-over regime. Since we are plotting the differ¬ 
ence the stabilizing effect of the far 

environment shows up as a growing positive hump. For 
each value of F, we plot three statistically independent 
numerical simulations. This gives us an idea about the 
statistical uncertainty of the results. We can clearly see 
that the curves which correspond to F = 0 are differ¬ 
ent from zero, due to the reasons discussed. Additional 
cases with increasing coupling to the heat bath. For those 
cases, the relative coupling strength F/A is 0.1 (circles), 
0.5 (triangles), and 1.0 (shaped crosses). We can observe 
that the theory agrees with the simulations, only in the 
case of smallest coupling, for stronger coupling the effect 
is systematically overestimated. In Fig. 2(b) we intend to 
find a rescaled decoherence rate Fgt, which best describes 
the numerical results, and hence the stabilizing effect of 
the far environment on the central system. A good agree¬ 
ment could be achieved only for times up to fmax ~ 15, 
which is the approximate location of the maxima of the 
curves shown. 

In Fig. 3 we repeat the comparison for A = 0.02, where 
the coupling between central system and RMT environ¬ 
ment is close to the perturbative regime. We use the same 
fitting procedure as in Fig. 2(b). Here, the values for the 
relative coupling strength F/A range from 0.1 to 5.0. The 
slowing down of decoherence in the central system due 
to the increasing coupling to the far environment, occurs 
just as in the case A = 0.1. However, for large values of 



FIG. 4. Fflt/A vs. F/A for A = 0.1, Ae = 25 and Ae = 50 
(crosses and squares respectively) and A = 0.02, Ae = 25 
(shaped crosses). The solid lines show the functions Fat/A = 
a and ba/{b + a) with b = 3.77 and a = V/\. 


F/A the deviations between simulations and theory be¬ 
come quite noticeable, even if we use the best fit values 
Ffit for the theory. 

Finally, we compare in Fig. 4 the fitted values Ffit for 
the coupling to the far environment, with the nominal 
ones, by plotting Fgt/A versus F/A. This is done for 
different dimensions of the near environment, for differ¬ 
ent coupling strengths between central system and near 
environment, and different couplings to the far environ¬ 
ment. The derivation of our theoretical result within 
linear response theory showed that the deviation from 
the exact result should be quadratic in F. The results 
for Ffit/A shown here, confirm this expectation, as they 
approach the line Ffit/A = F/A for sufficiently small val¬ 
ues. For larger values of F, the fitted values for Ffit, and 
thereby the stabilizing effect of the far environment, in¬ 
crease ever more slowly. To illustrate this behavior, we 
plotted the straight line Ffit/A = a as well as the func¬ 
tion g{a) = ba/{b+ a), with a = F/A and a best fit 
value b = 3.77, which describes the overall behavior of 
the points quite well. 

Summarizing, we have been able to obtain an ana¬ 
lytic expression confirming that nested environments can 
improve coherence of a central system as the coupling 
between near and far environment increases, as long as 
this coupling is small. We also extended previous limited 
numerical evidence for large coupling using a Caldeira- 
Leggett master equation which has been derived from 
RMT considerations in previous work [6]. This confirms 
that the effect subsists at large couplings between near 
and far environment, but subsides if the central system is 
strongly coupled to the near environment. An explana¬ 
tion on the basis of the quantum Zeno effect is tempting 
but problematic, at least in as far as we consider weak 
couplings between near and far environment. 
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Appendix: Derivation of the main result (Eq. 8) 

1. Dephasing coupling 

Under dephasing coupling, the nondiagonal element of the qubit reduced state is just the fidelity amplitude of 
the RMT-environment with respect to the perturbation induced by the coupling between central system and near 
environment. For an initial state ge,f, and with 14,f = UJ C) 14 and H\ = Hq -|- A14ff from Eq. (5) of the main article, 

= tr[£.e,f 

(A.l) 

The later two evolution operators are separable and therefore simplify as follows: 

= Mx{t) 0 If , Mx{t) = 

Since trf[AM 0 If] = Aij^ki Mkmhj = trf(A) M then 

= tre[ 0 e(i) Mx{t)], Qe{t) = tl'f [ ^-i(Ho+Ht)t/h^ 

(A.2) 


2. Linear response approximation for the coupling to the far environment 


The trace over the far environment in Eq. (A.2) is almost exactly of the form as the reduced density matrix (in the 
interaction picture) treated [11], namely with Q-i(Hx+Ht+'iVe^{)t/n^ we may write 

=tre[ge{t) Mx{t)] , Qeit) =tr{[Mx,'y{t) QeJ ■ (A.3) 

However, in order to apply the formalism of [11], we should assume the coupling 14,f and the initial state gie.f to be 
separable: 

l?e,f — He 0 ^ ^ 5 

and: 

h;^f(r) = dxir) 0 14(r) = Uj 0 Uf , ixir) = U,' , 

where we have already defined the representation 14,f(i) of the coupling operator to the far environment in the 
interaction picture and similarly for 14 (t). Of course there remains the very important difference, that here we 
have different echo operators on the left and the right side of the initial state. Nevertheless, following carefully the 
calculation in [11], developing the echo operators into their respective Dyson series we find: 

Qeit) = Qe - -^{Aj - Ai) , A j = Uf [ (t) He,f + He,f ^0 (i)^ ] , A/ = Uf [/a (t) He,f 4 (0^ ] 

>^A(i) =[ dr [ dr' 14 y(r) 14y(r') , hit) = f dr 14 ,f(r) . 


0 


0 


(A.4) 
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The calculation for the average over J\{t) with respect to the random matrix Vf yields 


= [ dr f dr' c{t - t') vx{t) vx{t') 'S)tf , 
Jo Jo 


(A.5) 


where c(r) describes the spectral correlations of Hf and (3 is the Dyson parameter, such that for a GUE (/3 = 2) or a 
GOE (/3 = 1) with Heisenberg time th = 27r?i/(io: c(t) = 3 — /3 + dir/i^) — &2(t/th). Similarly for Ai\ 



drdr' c(r 


t') vx{t) Qe Voir') . 


(A.6) 


Einally, we obtain 


{Aj - Al) = [ dr f dr' c{t - r) {vx{t) ['i)A(T0 ] - [ Va(t') ^o(t') ] 'i^o(T) } (A.7) 

Jo Jo 


3. Fermi golden rule regime and master equation 

If we assume that the Heisenberg time of the far environment th is very large, and that we are in the Fermi golden 
rule regime for the coupling to the far environment, then from the correlation function c(r), we only need to take the 
delta function into account. That reduces Eq. (A.7) to 

n t 

{Aj - Al) = ^ y dr ['DA(r) Vx{t) Pe - 2 Vx{t) Qe Vo{t) + Qe vo(t) vo(t) ] (A.8) 

Next, we will average that expression over the coupling matrix Ve, which is the near environment part of the coupling 
between near and far environment. Since this matrix is assumed to be an element of the GUE, we find: 

= X] ^ (^a(t) Vx(t)} = Ne (A.9) 

3 

On the other hand, we find 

Qe ^o{^))iq (^a)^ (.^e)jk {^x)kl Qlm (^e)np ('J'o)pq 

~ ('^x)ij ^kn 3jp {ux)kl Qlm (Uo)mn (uo)pq = {u\)ij {ux)kl Qlm {Uo)mk {uo)jq , (A.10) 

where ux = exp{—iHxt). This can be written as 

{vx{t) QeVo{T))^q = [Mx{T)^]iqtr[QeMx{T)] since Mx{t)=uIux. (A.11) 

Therefore, we obtain for ge{t) 

ge{t) = (^Net ge-2 J dr tr[pe A/a(t)] Ma(t)^^ (A.12) 

Let us denote E = 2tt Ng / (hdf), where we introduced the average level spacing df = h/rn- Then we obtain for the 
fidelity amplitude: 


/A.r(t) = tr (1 - rt) pe Mx{t) + — / dr tr[pe AfA(r)] Mxir)^ Mx(t) 


(A.13) 


thus 


/A.rW-^(l-rt)/A(t)+r f dr fxir) fx{t - t) (A.14) 

Jo 

where we have used that fx{t) = tr[pe Mx{t)] and Ng fx{t — T) = tr[MA(t — r)]. So fx{t) denotes the fidelity amplitude 
in the near environment, if there is no coupling to the far environment (j = 0). The first line is exact (in the limit 
rt <C 1), assuming that no ensemble averaging has been applied with respect to Jig and Vg. The second line assumes 
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self averaging for the quantities tr[pe Ma(t)] and tT[Mx{t — r)] which will probably hold for generic initial states Pe 
and sufficiently large near environment (Ng 3> 1). 
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